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We study the dynamics of Bose-Einstein condensates in tilted and driven optical superlattices. 
For a bichromatic lattice, each Bloch band split up into two minibands such that the dynamics is 
governed by the interplay of Bloch oscillations and transitions between the bands. Thus, bichromatic 
potentials provide an excellent model system for the study of nonlinear Landau- Zener tunneling and 
allow for a variety of applications in matter wave interferometry and quantum metrology. In the 
present paper we investigate the coherent dynamics of an interacting Bose-Einstein condensate as 
well as its stability. Different mechanisms of instability are discussed, which lead to a rapid depletion 
of the condensate. 

PACS numbers: 03.75.Lm,03.65.Sq 



I. INTRODUCTION 

Despite its apparent simplicity, the dynamics of quan- 
tum particles in periodic structures is full of surprises. 
Contrary to our intuition, a weak external field inhibits 
quantum transport in a periodic potential in favor of the 
celebrated Bloch oscillations |l|, |2j. For stronger fields 
a directed motion is re-introduced by repeated Landau- 
Zener transitions to higher Bloch bands [3|-[6| . One of the 
most interesting physical systems to explore the dynam- 
ics in periodic potentials are Bose-Einstein condensates 
(BEC) in optical lattices, allowing an in situ detection of 
the atoms @, Q . In recent years it became possible to 
realize periodic potentials with almost arbitrary shapes 
and an astonishing precision. Bichromatic lattices have 
been implemented by superimposing two incoherent op- 
tical lattices [H, [Toj| , or by combining optical potentials 
based on virtual two-photon and four-photon processes 
pj]-[T5|. These superlattices allow to engineer the Bloch 
band structure of the system by tuning few experimental 
parameters. 

The dynamics of a BEC in a tilted optical lattice, 
and especially the Landau-Zener tunneling between Boch 
bands, is strongly modified by the inter- atomic inter- 
actions. While the fundamental problem of Landau- 
Zener tunneling between two levels was solved as early 
as 1932 independently by Landau, Zener, Majorana and 
Stiickelberg [3|-l6| , a generalization of these results to in- 
teracting many-particle systems remains an open ques- 
tion up to today. The possibility to investigate Landau- 
Zener tunneling of a BEC in situ in a well controllable 
laboratory experiment has thus attracted much interest 
in recent years. A distinguished result of these studies 
was that strong interactions in a BEC lead to a break- 
down of adiabaticity and instability even in the limit of 
very slow parameter variations. This phenomenon was 
first predicted theoretically within a mean-field approxi- 
mation [l6l-[l9| and demonstrated for a BEC in an accel- 



erated optical lattice shortly afterwards [20l424j . Later, it 
was shown theoretically that the breakdown of adiabatic- 
ity results from the occurence of diabatic level crossings 
in the many-body spectrum [25l-[27|. 

In a bichromatic optical lattice, it is possible to tune 
the energy gap between the minibands via the relative 
phase of the two lattices and thus to control Landau- 
Zener tunneling as shown in [12|. Furthermore, the band 
stucture can be engineered such that the BEC is con- 
fined to the two lowest minibands and cannot be lost 
by repeated tunneling to higher excited bands, such that 
one can even observe repeated Landau-Zener tunneling 
events [15|, [28|, 29]. In addition, Landau-Zener tunneling 
can be used as a coherent beam splitter for atomic matter 
waves, which enables a variety of possible applications in 
matter wave interferometry and quantum metrology [HI . 

In the present paper we will provide a thorough the- 
oretical analysis of nonlinear Landau-Zener tunneling in 
bichromatic optical lattices with a focus on the onset of 
dynamical instability and depletion of the condensate. 
Furthermore, we study Bloch-Zener oscillations, the co- 
herent superposition of Bloch oscillations and Landau- 
Zener tunneling, and the dynamics in driven lattices, 
where transitions between the minibands are caused by 
the periodic driving instead of a static field. In the latter 
case, the dynamics of a BEC is governed by the interfer- 
ence of two matter waves with the same quasi momentum 
and opposite group velocity. This gives rise to an oscil- 
latory motion which can be viewed as a simple quantum 
simulation of the Zitterbewegung of a relativistic Dirac 
spinor. 

We focus on a bichromatic optical lattice with an al- 
ternating depth of the lattice wells. The many-body dy- 
namics in such an optical lattice is described by the Bose- 
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Hubbard type hamiltonian [3 1 
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where a n and are the bosonic annihilation and cre- 
ation operators, respectively. The parameter J denotes 
the tunneling rate between the wells, U is the interaction 
strength and F the strength of the static external field 
which accelerates the atoms. The parameter 5 > is the 
difference of the on-site energies between adjacent wells. 
It is directly proportional to the intensity of the double- 
periodic optical lattice. In the following we set J = 1 in 
all simulations, i.e. we measure all energies in units of 
the tunneling matrix element J. 

One of the main objectives of the present paper is 
a careful analysis of how the interactions affect the co- 
herent dynamics of a BEC in a bichromatic lattice and 
how they possibly lead to instabilities of the condensate. 
Throughout the paper we thus assume that the system 
is initially prepared as a pure BEC with N particles: 
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As long as the quantum state remains close to a pure 
BEC, the dynamics is well described within a mean- 
field approximation. The dynamics of the condensate 
wave function is then given by the discrete nonlinear 
Schrodinger equation (DNLSE) @ 



= -J(^n+l+^n-l) + 
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(3) 

However, dynamical instabilities lead to a depletion of 
the condensate such that the mean-field approximation is 
no longer applicable. In order to simulate the dynamics 
beyond mean-field we use the Bogoliubov backreaction 
(BBR) method which also gives a quantitative prediction 
for the depletion of the condensate (Hi- 



ll. BLOCH STATES AND BANDS 

Bloch states are the simultaneous eigenstates of the 
field-free Hamiltonian and a translation over two lattice 
sites. The dynamics of a BEC can be understood to a 
large extent from the properties of linear and non-linear 
Bloch states, including the depletion of the condensate. 
In the following we thus give a detailed analysis of the 
Bloch states and their stability. 
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FIG. 1: Nonlinear Bloch bands in a bichromatic optical lattice 
with 5 = 1 for g = -2 (a), g = +2 (b), g -0.5 (c) and 
g — +0.5 (d). Linear Bloch bands (g = 0) are shown in (e) 
for S = 1 (solid line) and S = (dashed line) for comparison. 



calculated as [28| 

E a {*)= V^ 2 +4J 2 cos 2 (^), (4) 

where a = 0, 1 labels the two minibands and k G 
[— 7r/2, +7r/2] denotes the quasimomentum. The band 
gap between the two minibands is directly given by the 
parameter 5. The corresponding Bloch states are given 
by \Xa iK ) = &LJ0), where 
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and 



u K = 4Jcos(ft) 
v K 



5+ V^+Ij^T cos 2 («), 
N K = 7r(uf. + v 2 ) being a normalization constant. 



A. Single-particle Bloch states 



B. Nonlinear Bloch bands 



In the simplest case of a single atom, where interac- 
tions are obviously irrelevant, the Bloch bands are easily 



Also in the mean-field approximation one can analyt- 
ically calculate the 'nonlinear' Bloch states, which are 
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defined as stationary states of the DNLSE 
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with the translation symmetry (/> n +2 = e 2m (/> n . Making 
the ansatz 
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one finds that the coefficients u K ,v K are determined by 
the two-mode DNLSE 
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Using the normalization \u K \ 2 + \v K \ 2 = 1, the effective 
coupling constant is given by g — 2U p, p being the aver- 
age particle density. Examples of nonlinear Bloch bands 
are shown in Fig. [TJ One observes that the bands be- 
come strongly asymmetric - for a repulsive nonlinearity 
g > the curvature of the ground band increases while 
the curvature of the excited band decrease and vice versa 
for an attractive nonlinearity g < 0. For strong nonlin- 
earities novel stationary states appear at the band edge 
k = ±7r/2, forming the so-called looped Bloch bands 
[34], [35|| . A quantitative analysis shows that DNLSE (jSJ) 
admits four solutions if [18| 
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Thus the critical nonlinearity for the existence of looped 
levels is directly linked to the band gap S. The defor- 
mation of the Bloch bands has significant consequences 
for the dynamics which will be discussed in detail in 
Sec. IIII Al An adiabatic dynamics is hindered by a sharp- 
ening of the levels and becomes completely impossible as 
soon as the loops form. 



C. Stability analysis 

The nonlinear Bloch states calculated in the preceding 
section are stationary states of the DNLSE (|3]). How- 
ever, they can become dynamically unstable due to the 
nonlinear interaction term, which also indicates a rapid 
depletion of the BEC. In order to determine the stability 
properties of a Bloch state flTJ), we add a small perturba- 
tion with quasimomentum q: 
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and analyze the consequences for the energy and the dy- 
namics of the condensate. If every perturbation increases 
the total energy given by the Gross-Pitaevskii energy 
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FIG. 2: Stability map for a nonlinear Bloch state with quasi- 
momentum k in the ground band for 5 = 1 and g — — 2 (a), 
g = +2 (b), g = -0.5 (c) and g = +0.5 (d). The gray scale 
map shows the stability exponent, i.e. the growth rate of a 
perturbation with quasimomentum q. Dynamical instability 
leads to a depletion of the condensate if the stability exponent 
is non-zero for at least one value of q. 
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then the Bloch state represents a local energy minimum 
and thus a stable superflow. Otherwise a perturbation 
may lower the energy and the Bloch state suffers a Lan- 
dau instability. In the present paper we are more con- 
cerned with the dynamical instability of a Bloch state, 
which occurs if a perturbation grows exponentially, as 
this indicates a rapid depletion of the condensate mode 
[36|-[39| . Note that dynamical instability always indicates 
energetical instability but not vice versa [35|. 

In order to determine the energetical stability we sub- 
stitute the ansatz ([TO]) into the Gross-Pitaevskii energy 
functional ([TT]) and expand it up to second order in the 
perturbation. The variation of the energy is then given 

by 
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Here we have introduced the abbreviations 

H (k\ = (I + 29 ^ 2 ~ » " 2Jcos ( fc ) "\ 
[ ) V -2Jcos(fc) -f + 2g\v\ 2 -fi] ' 

v= ( U o ° 2 ) and 

3« = (Co,K 5 Co, ^ 7 Ci,^) T - (14) 

The Bloch state represents a stable energy minimum if 
is positive for any perturbation, i.e. if the matrix 
L en (K,q) is positive definite for every q. 

The dynamical stability properties are found by substi- 
tuting the ansatz ([TO]) into the DNLSE (|3]). In first order, 
the perturbation evolves according to the Bogoliubov-de 
Gennes equation 

= L B dc(^, q)^K (15) 

with 
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FIG. 3: (Color online) Many-particle Landau-Zener tunneling 
between two minibands with band gap 5 = 0.2 (a) and 5 = 1 
(b). The dynamics has been calculated with the BBR method 
for the parameters J = 1, F = 0.2, UN 2 and N 100 
particles. The upper panels show the survival probability ([T5]) 
in the upper half of the lattice (solid line). 



L B dG(K,q) = cr z L en (n,q) 

A dynamical instability occurs if a perturbation grows 
exponentially, i.e. if there is any q for which the eigen- 
values of the matrix LedG^^) are n °t purely real. 

The dynamical stability for the Bloch states in the 
ground band is depicted in Fig. [2] for the same parame- 
ters as in Fig. [TJ A grey scale map shows the stability 
exponent, i.e. the maximum imaginary part of the eigen- 
values of the Bogoliubov-de Gennes matrix LBdG(ft, q) in 
dependence of ft and q. This imaginary part indicates 
the growth rate of a perturbation with wavenumber q 
and thus also the depletion rate of the condensate. A 
Bloch state with quasimomentum n is dynamically sta- 
ble only if the growth rate is zero for all values of q. One 
observes that two different kinds of dynamical instabil- 
ity exist in a bichromatic lattice. In the attractive case 
g < 0, the Bloch states in the center of the Brillouin 
zone become strongly unstable already for a quite mod- 
est nonlinearity. Thus we face the surprising conclusion 
that an attractive interaction on the one hand flattens 
the ground band and thus faciliates an adiabatic evolu- 
tion, but on the other hand leads to instability. On the 
contrary, a strong repulsive nonlinearity is required to 
introduce a weak dynamic instability at the edge of the 
Brillouin zone, which is associated with the occurence of 
looped Bloch bands. Thus one can infer that a significant 
depletion takes place only for a much stronger interaction 
than in the attractive case and that it sets in at the edge 
of the Brillouin zone around k = tt/2. 

III. DYNAMICS 

The dynamics of a BEC in a bichromatic optical lattice 
is governed by the interplay of the intraband dynamics 



and the transitions between the two minibands. A static 
force can be introduced by gravity [20|, magnetic field 
gradients J4l| , or by an acceleration of the complete lat- 
tice (see [8] and references therein). This force acceler- 
ated the atoms until they are finally Bragg reflected, lead- 
ing to the celebrated Bloch oscillations [TJ7 0, . In addi- 
tion, this force induces Landau-Zener tunneling between 
the minibands [12|, which will be investigated in detail 
in the following. In a bichromatic lattice one can control 
the tunneling rate to a large extent by tuning the lattice 
parameters [ll|. On longer timescale repeated Landau- 
Zener tunneling takes place, which leads to a complex 
dynamics due to the interference effects of atoms in the 
two minibands [15|, [28|, [29| . Finally, we will extend our 
analysis to the case where transitions are not driven by 
an external field but by a periodic driving. 

In order to simulate the dynamics beyond the mean- 
field approximation and to describe the depletion of the 
condensate we use the Bogoliubov backreaction (BBR) 
method introduced in [3ll-[33|] . In this approach one takes 
into account two- and four point functions and trun- 
cates all higher order correlation functions to obtain a 
closed set of evolution equations. The BBR method has 
proven its worth to predict the features of the many-body 
quantum state, especially the depletion of the condensate 
mode, avoiding the common problems of the Hartree- 
Fock-Bogoliubov approximation [3l|-[33j. In particular 
the nature of the many-body quantum state is indicated 
by the reduced single-particle density matrix (SPDM) 

<*j,m = -^(^m)- (17) 

The leading eigenvalue of this matrix, Aq, gives the frac- 
tion of atoms in the condensate mode. Consequently, the 
non-condensed fraction is given by 1 — Aq. 
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band gap 5 

FIG. 4: (Color online) Landau-Zener tunneling probability 
Plz between two minibands as a function of the band gap 5 
for UN = -2 (o), UN = (•), UN = +2 (o) and UN = 
+4 (v) and N = 100 particles. The solid line shows the 
analytic approximation (|19|) for the linear case. Dashed lines 
are plotted to guide the eye. 



A. Nonlinear Zener tunneling 

To begin with, we explore the basic features of the 
dynamics in a bichromatic lattice subject to a static ex- 
ternal field, Bloch oscillations and Landau-Zener tunnel- 
ing, for the weakly interacting case. We assume that 
the initial state is a pure BEC (j2j) in the ground Bloch 
band with k = weighted by a Gaussian envelope 
ip n (0) ~ (/> n exp(— (n — no) 2 /4<j 2 ) with a width of a = 5 
sites centered around the site no = 35. In the forth- 
coming examples we choose the total particle number to 
be N — 100 located in a lattice with M = 46 sites and 
F = 0.2, unless otherwise stated. A weak nonlinearity 
induces a reversible dephasing, which damps Bloch oscil- 
lations (39I44H. and, of course, alters the Landau-Zener 
tunneling rate between the two minibands. 

Figure [3] shows two examples of the many-body dy- 
namics starting from the initial state (j2j) for a weak re- 
pulsive interaction, UN = 2, and two different values of 
the band gap 5. The figures show the evolution of the 
density (hj(t)), j = 1,...,M in false color. Here and 
in the following, we take the Bloch time Tb = Ztt/F of 
the single-periodic lattice (5 = 0) as the reference time 
scale. One observes that the BEC is first accelerated by 
the external field F until it reaches the edge of the Bril- 
louin zone at T = Tg/4. If the band gap 6 is large (Fig. [3] 
(b)), the BEC matter wave stays in the ground miniband 
and performs Bloch oscillations with a period of Tg/2. 
In contrast, if the gap is small, the matter wave tunnels 
to the excited miniband and performs Bloch oscillations 
with the full period T B (Fig. El (a)). At time t = T B /2 it 
is located at the turning point of the Bloch oscillations. 
For intermediate values of the band gap, only a fraction 
of the condensate tunnels to the excited miniband and 
the wavepacket splits. 

For a further quantitative analysis of the Landau-Zener 
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FIG. 5: (Color online) Unstable dynamics in a tilted bichro- 
matic optical lattice for strong interactions UN = +20 (a,b) 
and UN = —10 (c,d). The left panels (a,c) show the evolution 
of the density and the right-panels (b,d) show the magnitude 
of the scaled SPDM N\a 3 , k \ = \(a]a k )\ at t = T B /2 in a 
colormap plot. The remaining parameters are 8 = 0.2 and 
N 1000. 

tunneling rate we estimate the survival probability in the 
ground miniband by the number of atoms remaining in 
the upper half of the lattice: 

1 46 

P svW = ^E<^)>' ( 18 ) 

The Landau-Zener tunneling probability to the excited 
band is then given by Plz := 1 — Psv(Tb/2). The time 
dependence of the survival probability ([18]) is shown in 
the upper panels of Fig. [3] together with Plz- 

Figured] shows the Landau-Zener tunneling probability 
Plz in dependence of the band gap S for different values 
of the interaction strength UN. In the linear case, UN = 
0, one can approximate the avoided crossing of the two 
minibands at the edge of the Brillouin zone by an effective 
two-level model which yields the following approximation 
for the Landau-Zener probability [28|: 

This approximation shows an excellent agreement with 
the numerical results shown in Fig. 2J In the weakly non- 
linear case one observes an increase of the Landau-Zener 
tunneling rate Plz for a repulsive nonlinearity U > and 
a decrease for an attractive nonlinearity U < 0, which has 
also been demonstrated experimentally [2l|. This effect 
can be understood from the structure of the nonlinear 
Bloch states introduced above. With increasing inter- 
action strength, the nonlinear Bloch bands become 
strongly asymmetric as shown in Fig.[TJ For UN < 0, the 
ground band is flattened so that adiabaticity is faciliated 
and Plz decreases, while the excited band is sharpened. 
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FIG. 6: (Color online) Depletion of the condensate during 
nonlinear Landau- Zener tunneling. Temporal growth of the 
non-condensed fraction 1— Ao for a strong repulsive interaction 
UN — +20 (a,b) and an attractive interaction UN = — 10 
(c,d). The remaining parameters are 6 = 0.2 (a,c) and S = 1.2 
(b,d), respectively, and N = 1000. 




FIG. 7: (Color online) The non-condensed fraction 1 — Ao at 
time t — Tb/2 as a function of (a) the interaction strength 
UN for N = 100 (o), N = 300, (□) and N = 1000 (o) and (b) 
as a function of the particle number N for a fixed interaction 
strength UN = +6 (o) and UN = -6 (□). Solid lines are 
drawn to guide the eye. 



The inverse effect is found for UN > such that Plz in- 
creases. If the nonlinear ity UN exceeds a critical value, 
a Bloch state at the edge of a band bifurcates to a looped 
structure, which prevents an adiabatic evolution even for 
very small values of F. 

B. Depletion of the condensate 

Strong inter-atomic interactions alter the dynamics of 
the BEC completely. Examples are shown in Fig. [5] for 
a repulsive (a,b) and an attractive (c,d) interaction, re- 
spectively. One observes that the familiar Bloch oscilla- 
tion pattern is significantly disturbed, especially in the 
case of attractive interactions. In the repulsive case, the 
atoms are distributed over several lattice sites, but the 
phase coherence between theses sites is lost almost com- 
pletely. This is indicated by a strong suppression of the 
non-diagonal parts of the SPDM as shown in Fig. [5] (b). 
A strong attractive nonlinearity leads to a collapse of 
the condensate. Figure [5] (c) shows that the atoms are 
strongly focussed to a single lattice site at t « 0.15 Tg. 
Afterwards, a fraction of the atoms 'explodes' from the 
focus and the condensate mode is rapidly depleted. 

To further analyze the different mechanisms of instabil- 
ity due to repulsive and attractive interactions we calcu- 
late how the condensate is depleted. Figure [6] shows the 
time evolution of the non-condensed fraction 1 — Ao for 
a strong repulsive (UN = +20) and a strong attractive 
nonlinearity (UN = —10), respectively. In the attractive 
case, instability sets in much earlier and takes place on 
a very short time scale. This difference can be explained 
by the results of a linear stability analysis as discussed 
in Sec. Ill CI In the repulsive case the onset of dynam- 
ical instability can be associated with the emergence of 



looped Bloch bands. The condensate becomes dynam- 
ically unstable at the edge of the Brillouin zone where 
the loops emerge. In contrast, already modest attractive 
interactions lead to a dynamic instability at the center of 
the Brillouin zone (cf. Fig [2j) such that the depletion of 
the condensate sets in immediately. 

A quantitative analysis of the depletion of the conden- 
sate is provided in Fig. [7J where we have plotted the 
non-condensed fraction at t = Tb/2 as a function of the 
interaction strength UN in (a) and the particle number 
N for a fixed value of the interaction strength UN = ±6 
in (b). Figure 0(a) clearly shows the qualitative differ- 
ence between an attractive and a repulsive interaction. In 
the first case, one observes a rapid increase of the non- 
condensed fraction when the interaction strength exceeds 
the critical value for the onset of a dynamical instabil- 
ity. For a repulsive interaction, however, the dynamics 
is rather stable so that the non-condensed fraction re- 
mains small for all values of \UN\ < 10 shown in the 
figure. The non-condensed fraction decreases with the 
particle number and tends to zero in the mean-field limit 
N — )> oo. However, the speed of convergence depends cru- 
cially on the stability of the dynamics as shown in Fig. 
(b). In the repulsive case, UN = +6, the dynamics is 
stable and thus convergence is fast. The non-condensed 
fraction decreases rapidly with increasing particle num- 
ber such that the mean-field description by the DNLSE 
is valid already for quite small values of the particle num- 
ber. On the contrary, the convergence is logarithmically 
slow for UN = —6 due to the dynamical instability. A 
different approach to a classical instability causing de- 
pletion is provided by generalized mean-field descriptions 
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C. Bloch-Zener oscillations 

On a longer timescale, the dynamics of a BEC in a 
tilted optical lattice is governed by the interference of 
Bloch oscillations and Zener tunneling between the Bloch 
bands [28|, • Figure [8] (b) shows an example of the 
dynamics of the atomic density for S = 0.5. The con- 
densate wave packet is coherently split by Landau- Zener 
tunneling between the two minibands at t = Tg/4 and 
recombined again at t = 31#/4, thus realizing an ef- 
fective matter wave Mach-Zehnder interferometer. The 
splitting ratio of this interferometer, which is given by 
the Landau- Zener tunneling rate (fT9]h is easily tunable 
by changing the band gap S. 

For very small and for very large values of the con- 
densate occupies only one miniband - it remains in the 
ground band for large S and tunnels completely to the 
other miniband for small S as shown in Fig. [8] (c). In 
both cases the condensate shows simple Bloch oscilla- 
tions and returns back to its initial state at t = Tb- For 
intermediate values of 5, the condensate is split into two 
parts at t = Tb/^- The further dynamics and especially 
the occupation of the two minibands is governed by the 
interference of the two possible paths. For the given pa- 
rameters, about one half of the population is still local- 
ized in the excited miniband at t = Tb = 2tt/F. In this 
parameter range, the dynamics is very sensitive even to 
small nonlinearities as shown in Fig. [8](d). The survival 
probability at t = Tb differs significantly for UN = —1 
and UN = +1, although the nonlinearity is still compar- 
atively weak. 

For the given interaction strength \U\ < 1, the splitting 
and also the recombination of the condensate is fully co- 
herent; the fraction of non-condensed atoms is less than 
0.4% at t = 2Tb as shown in Fig. [8] (a). A significant 
depletion of the condensate is observed only for stronger 
nonlinearities; for instance the non-condensed fraction at 
t = 2T B exceeds 10% for UN > 5. 

This example demonstrates the possible use of Landau- 
Zener tunneling and Bloch-Zener oscillations in quantum 
metrology. These tools can be used, for instance, to di- 
rectly measure the band structure of a bichromatic po- 
tential as demonstrated in [15]. This is a unique feature 
of bichromatic optical lattices. In a simple periodic po- 
tential, a matter wave will be accelerated further towards 
— oo after it has escaped from the ground band, such that 
no interference can be observed. 



D. Coupling of bands by a periodic driving 

Previously, we have discussed the effects of Zener tun- 
neling between the two minibands induced by the ex- 
ternal field F. A coupling of the bands can also be 
introduced in the field free case by a periodic driving 
of the system parameters. This has the advantage that 
the quasi momentum n is conserved such that a different 
regime of the dynamics can be explored. 
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FIG. 8: (Color online) Bloch-Zener oscillations of a BEC with 
N = 300 atoms in a tilted bichromatic lattice, (a) Non- 
condensed fraction 1 — Ao for 5 — 0.5 and UN — 1. (b) 
Atomic density in real space for 5 — 0.5 and UN = 1. (c,d) 
Survival probability (|18f) in the upper half of the lattice after 
t = T B /2 and t = T B , respectively, for UN = -1 (o) and 
UN = +1 (o). 



Here we consider a BEC initially prepared in the 
ground miniband with a well defined quasimomentum k. 
The strength of the double-periodic optical lattice is var- 
ied in time to realize a harmonic driving of the energy 
offset 



S(t) = Sq + Si cos(u;£). 



(20) 



This driving induces transitions between the two mini- 
bands if the frequency is chosen to be resonant with the 
band gap, uj = E\(k, Sq) — Eq(k,8q). In the following 
example we set So = 0.4, S\ = 0.2 and U = 0. The 
initial state is assumed to be pure BEC with momen- 
tum n = 0.17T, weighted by a Gaussian envelope with 
(7 = 10. The resulting dynamics is shown in Fig. [9] in 
real (left) and momentum space (right). One clearly ob- 
serves the transitions between the two minibands, while 
the quasimomentum of the BEC is conserved (panel (b)). 
A further quantitative analysis of this effect is provided 
in panel (d), where the occupation of the two minibands 
Po y i is plotted. The oscillation between the bands has 
remarkable consequences for the real-space dynamics of 
the BEC shown on the left-hand side of the figure. As the 
two minibands have opposite curvature, a transition be- 
tween the bands reverses the group velocity of the matter 
wave. This leads to an oscillatory motion of the mean po- 
sition, which can be understood as a quantum simulation 
of the Zitterbewegung of a Dirac spinor. This relativistic 
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FIG. 9: (Color online) Quantum dynamics of BEC in a Dichro- 
matic lattice with periodic driving 8(t) and F = 0. (a,b) 
Atomic density in real and momentum space, respectively, 
(c) Position expectation value . j{ftj) /N. (b) Occupation 
of the two minibands po (blue) and pi (red). 

effect results from the interference of particle and anti- 
particle contributions moving to opposite directions. In 
the discussed quantum simulator, the two minibands thus 
play the role of particle and anti-particle contributions, 
respectively. Similar effects were recently predicted for 
optical waveguide arrays (42|, |43| . 

IV. CONCLUSION AND OUTLOOK 

We have discussed the dynamics of a Bose-Einstein 
condensate in a bichromatic optical lattice. In such a 
lattice, the ground Bloch band splits up into two mini- 
bands with a controllable band gap. Hence, they are 
ideally suited to study the complex quantum dynamics 
resulting from the interplay of the intraband dynamics 
and transition between the minibands. 

The basic dynamics of a BEC in an optical lattice 
can be inferred from the band structure of the system. 
Within the tight-binding approximation, one can readily 
calculate the linear as well as the nonlinear Bloch states. 
In particular, this treatment yields an explicit expression 
for the critical interaction strength for the occurrence 
of looped Bloch bands, which leads to a breakdown of 



adiabatic motion. The stability properties of the Bloch 
states have been analyzed in detail by the Bogoliubov-de 
Gennes approach. 

The dynamics of a BEC was simulated using the 
Bogoliubov backreaction method which also provides a 
quantitative estimate for the depletion of the conden- 
sate. The Landau-Zener tunneling of a BEC between 
minibands in a tilted or accelerated bichromatic lattice 
has been intestigated in detail. For weak interactions, 
the condensate remains essentially intact, whereas the 
Landau-Zener tunneling process is strongly affected. Re- 
pulsive interactions increase the tunneling rate in partic- 
ular in the 'adiabatic regime' of large band gaps, while at- 
tractive interactions suppress Zener tunneling. Strong in- 
teractions cause a dynamical instability and thus a rapid 
depletion of the condensate mode. However, the mech- 
anism of dynamical instability is significantly different: 
In the repulsive case, instability sets in at the edge of 
the Brillouin zone and is intimately related to the occur- 
rence of looped Bloch bands. A condensate with attrac- 
tive interactions is unstable already in the center of the 
Brillouin zone, leading to a collapse and explosion of the 
condensate. 

On longer timescales, the interplay of Bloch oscillations 
and Landau-Zener tunneling leads to a complex dynam- 
ics due to the interference of the contributions in the two 
minibands. A quantitative analysis of this effecs has been 
given, which also shows the applicability of Bloch-Zener 
oscillations in matter wave interferometry. A remark- 
able dynamics is also observed if transitions between the 
bands are not induced by a static external field, but by 
a periodic driving which leaves the quasimomentum un- 
changed. Because of the different curvature of the mini- 
bands, the matter waves in the two minibands move into 
opposite directions. The interference of the two contri- 
butions then leads to a dynamics which is comparable to 
the Zitterbewegung of a Dirac spinor. 
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